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Abstract 

We study the collision of two flat, parallel end-of-the- world branes in heterotic 
M-theory. By insisting that there is no divergence in the Riemann curvature as 
the collision approaches, we are able to single out a unique solution possessing the 
local geometry of (2d compactified Milne)/Z2 x M3, times a finite-volume Calabi- 
Yau manifold in the vicinity of the collision. At a finite time before and after the 
collision, a second type of singularity appears momentarily on the negative-tension 
brane, representing its bouncing off a zero of the bulk warp factor. We find this 
singularity to be remarkably mild and easily regularised. The various different 
cosmological solutions to heterotic M-theory previously found by other authors are 
shown to merely represent different portions of a unique fiat cosmological solution 
to heterotic M-theory. 
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1 Introduction 



The riddle of the initial singularity is one of the most basic challenges in cosmology. In 
standard four- dimensional general relativity, the Riemann curvature diverges at the big 
bang signalling an irretrievable breakdown of the theory. In higher-dimensional string and 
M-theory, however, the nature of the initial singularity is significantly altered. Within 
the higher-dimensional picture, the Riemann curvature may remain bounded all the way 
to the singularity. In this situation, string and M-theory corrections to the background 
geometry remain small, allowing one to attempt to study the propagation of strings and 
branes right up to, and perhaps even across, the singularity. 

In this paper, we wish to study such a model of the cosmological background spacetime, 
within an especially well-motivated theoretical framework. Heterotic M-theory is built on 
the profound correspondence between supergravity and strongly-coupled heterotic string 
theory [1, 2], and it remains one of the most promising approaches to the unification 
of particle physics and gravitation. From the perspective of particle phenomenology, 
heterotic M-theory makes full use of branes and of the extra M-theory dimension, in 
order to solve the fundamental puzzles of chirality and of the difference between the GUT 
and Planck scales. From a cosmological perspective, heterotic M-theory provides a setup 
in which standard model matter is localised on branes, allowing the background density 
of such matter to remain finite at the initial singularity. 

We consider the coUision of two flat, parallel end-of-the- world branes within heterotic 
M-theory [3, 4]. Prom the standpoint of the four-dimensional effective theory, this event is 
a cosmological singularity of the usual catastrophic type. Yet from the higher-dimensional 
standpoint, the situation is far less singular. The density of standard model matter 
remains finite. Moreover, even the bulk spacetime between the branes is relatively well- 
behaved: as we shall show, the Riemann curvature remains finite for all times away from 
the collision event itself. In the case of trivial, toroidal compactifications of M-theory, this 
phenomenon is well known. Here, the coUiding-brane spacetime is locally fiat, a product of 
two-dimensional compactified Milne space-time, modded out by Z2, with nine-dimensional 
fiat space. At the collision itself, the spacetime is non-Hausdorff since one dimension 
momentarily disappears. Nevertheless, the low-energy degrees of freedom, described near 
the collision by winding M2 branes (or, equivalently, weakly coupled heterotic strings), 
possess a regular evolution across the singularity [5]. 

In more realistic models, six spatial dimensions are compactified on a Calabi-Yau 
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manifold. The dynamics are then described by a five- dimensional efi^ective field theory 
known as heterotic M-theory, which is a consistent truncation of eleven-dimensional su- 
pergravity coupled to the boundary branes [6, 7]. Recently, an improved formulation has 
been developed which avoids problematic terms involving squares of delta functions [8, 9], 
giving one greater confidence that classical solutions of the five-dimensional effective the- 
ory do indeed provide consistent M-theory backgrounds. In this paper, we shall show 
that these equations possess a unique global solution representing the collision of two fiat 
end-of-the-world branes which, in the vicinity of the collision, reduces to (2d compactified 
Milne)/Z2 x R3, times a finite- volume Calabi-Yau manifold. In this solution, the Riemann 
curvature is again bounded at all times away from the collision event itself. Our solution 
offers the intriguing prospect of modelling the big bang as a brane collision in a setup 
with a high degree of physical realism. 

The idea that the big bang was a brane collision in heterotic M-theory was first 
proposed in [3, 4]. However, the solution representing the approach and collision of 
two branes has not so far been given in complete form. Important steps towards such a 
solution were taken by Chamblin and Reall [10], who found a static solution for the bulk 
geometry in which moving, spatially flat branes can be consistently embedded. The bulk 
geometry possesses a timelike virtual naked singularity lying 'beyond' the negative-tension 
boundary brane. Two qualitatively different solutions then exist according to whether the 
boundary branes move through the static bulk in the same, or in opposing, directions. 
Both of these solutions, together with their time-reverses, are naturally combined within 
our global solution. In a more recent paper [11], Chen et al. considered the solution 
where both the positive- and negative-tension boundary branes move towards the virtual 
singularity. They found an exact solution in a convenient alternative coordinate system 
which is comoving with the branes. However, Chen et al. found that both branes encounter 
the naked singularity, resulting in the annihilation of the entire spacetime. 

In this paper, we shall argue for a different fate of the Chen et al. solution. We 
shall show that the naked singularity is 'repulsive' as far as the negative-tension brane 
is concerned, so the brane approaches it with vanishing speed. The ensuing singularity 
is extremely mild and easily regulated, for example by introducing an arbitrarily small 
amount of matter on the negative-tension brane. The result is a smooth bounce, at flnite 
Calabi-Yau volume, so that the negative-tension brane recoils from the singularity and 
collides with the positive-tension brane. The end-of-the-world brane collision is locally 
(2d compactified Milne)/Z2 x M3, with a finite-volume Calabi-Yau manifold. Assuming 
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that M-theory can deal with this second singularity in the manner described by [5] , we can 
follow the system through. The M-theory dimension then re-appears and the negative- 
tension branc is thrown back towards the naked singularity. A second bounce of the 
negative-tension brane occurs before the system continues into the future with both branes 
expanding. The complete solution is illustrated in the Kruskal diagram shown in Figure 1. 

We will present the single, global solution in two different coordinate systems, each of 
which has some merit. In the first system, which comoves with the branes so that they are 
kept at fixed coordinate locations, only the bulk is dynamic. Here, we are able to derive 
the solution as a series expansion in the relative rapidity of the branes at the collision. 
This method has been previously apphed to solve for the bulk geometry and cosmological 
perturbations of the Randall-Sundrum model [12] . The chief virtue of this brane-comoving 
coordinate system is that it simplifies the junction conditions on the branes, and allows 
for a manageable treatment of cosmological perturbations. 

In the second coordinate system, the branes move whereas the bulk is static. The 
boundary conditions that we impose at the moment of the collision, in conjunction with 
the assumed cosmological symmetry and spatial flatness of the branes, allows us to derive 
a modifled Birkhoff theorem singling out a unique solution. Away from the brane collision, 
this solution has finite Riemann curvature throughout, at least after suitable regularisation 
of the bounce of the negative-tension brane. The ensuing evolution, before and after the 
latter event, is described by Chamblin and Reall's two solutions [10], the second of which 
Chen et al. succeeded in re-writing in a brane-comoving coordinate system. 

The outline of this paper is as follows. In Section 2, we review the standard static 
solution of heterotic M-theory and fix notation. In Section 3, we derive the colliding 
brane solution in brane-comoving coordinates, as an expansion in the rapidity of the 
brane collision. In Section 4, we discuss the same solution in a coordinate system in 
which the bulk is static and the branes are moving. Considering the relevant modified 
Priedmann equation, we then show how the negative-tension brane bounces when an 
arbitrarily small amount of matter is present on the branc. Three appendices are devoted 
to technical matters. The first proves a Birkhoff-like theorem for the bulk. The second 
derives the equations governing the motion of branes in the bulk, and the third discusses 
the Kruskal extension of the bulk geometry. These calculations are used to make accurate 
plots of the brane trajectories, like that shown in Figure 1. 
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Figure 1: A Kruskal plot of the entire solution, in the spirit of Chen et al. [11]. The exact 
trajectories of the positive- and negative-tension branes are plotted in red and green respectively, 
while the naked singularity is indicated by a thick black line. Representative orbits of the bulk 
Killing vector field are shown with dashed lines, while the solid straight lines indicate the Boyer 
axes. The bounces of the negative-tension brane off the naked singularity, as well as the collision 
of the branes themselves, are shown at a magnified scale in the inset. In this plot, we have chosen 
the relative rapidity of the brane collision to be 2yo = 1- Analogous plots for greater collision 
rapidities may be found in Figure 4 of Appendix C. 
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2 Heterotic M-theory 



Eleven-dimensional supcrgravity can be compactified on a Calabi-Yau 3-fold to give a min- 
imal five-dimensional supergravity theory [1, 2, 6, 7]. Although the dimensional reduction 
of the graviton and the 4-form flux generates a large number of fields, it is consistent to 
retain only the five- dimensional graviton, and a scalar parameterising the volume of 
the Calabi-Yau manifold (namely Vcy = e"^). This theory can be consistently coupled 
to two four-dimensional boundary branc actions, provided that one keeps, in addition 
to gravity and the scalar, the components of the 4-form flux with all indices pointing in 
the Calabi-Yau directions. This flux thus appears as a scalar in flvc dimensions (albeit a 
constant one, as can be seen from its Bianchi identity [6, 7]), and it leads to a potential 
for (j). The resulting dimensionally reduced effective action is given by 

Jbd ^ 



+12a 




where a is related to the number of units of 4-form flux^, and where we have placed 
branes of opposite tensions at y = ±yo {y being the coordinate transverse to the branes) . 
In brane-comoving coordinates, the resulting equations of motion are 

Gab = ^0,a0,& - ^gah(t>'''(t>,c - ^Qab 0?e~'^'^ 

+6a[<5(y + i/o) - b{y - y^)\b^Xg^,^ e't (2.2) 

y9yy 

□0 = -12a''e-"f' + 12a[5{y + yo)-5iy-yo)]^e-'(', (2.3) 

\/9yy 

where Latin indices run over all flve dimensions, but Greek indices run over only the four 
dimensions to which the y coordinate is normal. The static domain wall vacuum solution 
is given by 

ds^ = h^/\y){7]^,dx^dx'' + dy^), 

h{y) = ba{y + yo) + c, -yo < y < yo- (2.4) 

The y coordinate is taken to span the orbifold S^/Z2, with fixed points at ±yo. In 
an 'extended' picture of the solution, obtained by Z2-reflecting the solution across the 

^Compared to [6, 7], we have rescaled a such that a = (aLOSw/Sv^). 
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branes, there is a downward-pointing kink at y = —yo and an upward-pointing kink at 
y = +1/0 • These ensure that the appropriate Israel matching conditions are satisfied, with 
the negative-tension brane being located at y — — yo and the positive-tension brane at 
y — +yo- The integration constant c is arbitrary. 

3 A cosmological solution with colliding branes 

We now wish to generahse the static solution above to allow for time dependence. In 
particular, we are interested in cosmological solutions in which the two boundary branes 
collide, with the five-dimensional spacetime geometry about the collision reducing to (2d 
compactified Milne) /Z2 x R^. We will therefore take as our metric ansatz 

ds^ = n\t,y){-df + t^dy^) +b\t,y)dx^, (3.1) 

where the span the three-dimensional spatial worldvolume of the branes, which we are 
assuming to be fiat. This ansatz is in fact the most general form consistent with three- 
dimensional spatial homogeneity and isotropy, after we have made use of our freedom 
to write the (t, y) part of the metric in conformally fiat Milne form. In the above, and 
throughout this section, we will take the branes to be fixed at the constant coordinate 
locations y — ±yo- 

The functions n and 6, as well as the scalar field 0, are arbitrary functions^ of t and 
y. In order that the collision be the 'least singular' possible, we require that the brane 
scale factors, as well as the Calabi-Yau volume, be finite and non-zero at the coUision. 
By an appropriate choice of units, we will therefore demand that n, h and all tend 
to unity as t ^ 0. In this way the geometry about the collision reduces to the desired 
(2d compactified Milne)/Z2 x R'^ times a finite-volume Calabi-Yau manifold. Since the 
five-dimensional part of this geometry is flat, our solution will be protected from higher- 
derivative string or M-theory corrections in the vicinity of the collision. 

To solve the relevant equations of motion, (2.2) and (2.3), we will perform a perturba- 
tive expansion in the rapidity of the brane collision. This latter quantity is given by 2yo 
(see e.g. [14]), and so we are interested in the case in which yo ^ 1? ^-c, the case in which 
the relative velocity of the branes at the collision is small. (See [12] for an analogous 
procedure in the case of the Randall-Sundrum model). 

^Note in particular that we are not assuming a factorised metric ansatz {e.g. n{t,y) = ni{t)n2{y)), 
since this generically leads to brane collisions in which the Calabi-Yau manifold shrinks to zero size at 
the moment of collision [13]. 
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To implement this expansion, we first introduce the re-scaled time and orbifold coor- 
dinates 

uj = y/yo, T = y^t, (3.2) 

so that the branes are now located at a; = ±1, independent of the rapidity of the brane 
collision. For convenience, we will also adopt the convention that primes denote derivatives 
with respect to u and dots denote derivatives with respect to InT, i.e., 

' = d^, ■ = Tdr. (3.3) 

In order to write down the junction conditions and the Einstein equations, it is useful to 
work with the variables = and e'^ — yont — nT. The junction conditions [15], vahd 
when evaluated ata; = —1"*" ora; = -|-l~, are then 

u'^ae"-'^, (5' ^2,a (j)' ^ 6a e"-'^ . (3.4) 

Evaluating the 0, Gq + G^, G^, G^ + G^ — (1/2)G-, and Go5 equations in the bulk, we 
obtain the following set of equations: 

0" + + 12a V^-^-^ = yl{4> + /3<j)), (3.5) 
P^ + P'^ + Qa^e""-'^ = yUP + h, (3-6) 
l/3'2^/3V-l</,'2 + 3aV^-^^ = yl0+lp^-Pu+\4>'), (3.7) 

+ = y'oii^-l$' + -J'), (3.8) 

$' + ^P(3'-0(3'-u'$+^,j><l>' = 0. (3.9) 

In the above, both the G^ equation (3.7) and the Gqs equation (3.9) involve only 
single derivatives with respect to uj. Applying the junction conditions, we find that both 
left-hand sides vanish when evaluated on the branes. The Gqs equation is then trivially 
satisfied, while the G| equation yields the relation 

valid on both branes. Introducing the brane conformal time r, defined on either brane 
via byodr = ndT, this relation can be re-expressed as 

Y + i^-O. (3.11) 
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Notice that in the Einstein equations (3.5)-(3.8), all terms involving time derivatives 
appear at 0{iIq) higher than the terms involving derivatives along the orbifold direction. 
Thus, solving the equations of motion perturbativcly in powers of y^, at zcroth order we 
have a set of ordinary differential equations in uj. Integrating these differential equations, 
we obtain the uj dependence of the solution at leading order, along with three arbitrary 
functions of time. How these functions are determined will be explained below in Section 
3.2. Then, at next-to-leading order, we have a set of ordinary differential equations for 
the UJ dependence of the solution at O(yo)) involving source terms constructed from the 
time dependence at leading order. Thus, solving the Einstein equations is reduced to 
an iterative procedure involving the solution of a finite number of ordinary differential 
equations at each new order in y^. 

3.1 Initial conditions 

In order to fix the initial conditions for our expansion in i/q, it is useful to solve for the bulk 
geometry about the collision as a series expansion in t (in this section we revert briefly to 
our original t and y coordinates). Up to terms of order i^, the solution corresponding to 
the Kaluza-Klein zero mode is: 

n — a (sech|/o sinh y)t-\ sech^|/o (9 — cosh 2yo — 8 cosh 2y) t^, (3.12) 

8 

b — (sechyo sinh y)t^- — sech^yo(3 + cosh 2yo — 4 cosh 2y) t^, (3.13) 

3q;2 

— 1 + &a (sechyo sinh y)t — sech^7/o(2 — cosh 27/o — cosh 2y) t^, (3-14) 

where we have used the junction conditions to fix the arbitrary constants arising in the 
integration of the bulk equations with respect to y. The brane conformal times t± (where 
the subscript ± refers to the brane locations uj — ±1, or equivalently, to their tensions) 
are then given by 

T± = y ^dt = t + 0(t'), (3.15) 

in terms of which the brane scale factors h± arc 

3 

h± = 1 ± Q;r±tanhyo - - Q;V| tanh^ yo + 0(t±)^ (3.16) 

3.2 A conserved quantity 

As explained above, every time wc integrate the bulk Einstein equations at a given order 
in t/q with respect to uj we pick up three arbitrary functions of time. Two of these three 
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functions may be determined with the help of the G| equation evaluated on both branes, 
namely (3.11). To fix the third arbitrary function, however, a further equation is needed, 
which we derive below. 

Introducing the variable x = — 2/3, upon subtracting twice (3.6) from (3.5) we find 



which is the massless wave equation = in this background.^ Since the junction 
conditions imply x' = on the branes, the left-hand side vanishes upon integrating over 
uj. A second integration over InT then yields 



for some constant 7, which we can set to zero since our initial conditions are such that 
Xe^ ^ as T ^ 0. 

Let us now consider solving (3.17) for x, as a perturbation expansion in y^. Setting x = 
Xo + 1/0 Xi + 0(yo)' similarly for at zeroth order we have (Xo^'^")' ~ 0- Integrating 
with respect to uj introduces an arbitrary function of T which we can immediately set to 
zero using the boundary condition on the branes, which, when evaluated to this order, 

read Xo = 0- This tells us that Xo = throughout the bulk; Xo is then a function of T 
only, and can be taken outside the integral in (3.18). Since 7 = 0, yet the integral of 
across the bulk cannot vanish, it follows that xo must be a constant. 

At order i/q, the right-hand side of (3.17) evaluates to yo(Xoe*)'; which vanishes. 
Evaluating the left-hand side, we have (x'lC^")' = 0, and hence, by a sequence of steps 
analogous to those above, we find that Xi must also be constant. It is easy to see that 
this behavior continues to all orders in y^. We therefore deduce that x = — 2/3 is exactly 
constant. Since both and /? tend to zero as T ^ 0, this constant must be zero, and so 
we find 



The essence of this result is that a perturbative solution in powers of y1 exists only 
when X is in the Kaluza-Klein zero mode. This may be seen from (3.17), which reduces, 
in the limit where T and /i^ ^ 0, to x" = VqX- "^^^ existence of a perturbative 
expansion in requires the right-hand side of this equation to vanish at leading order. 

•^For branes with non-zero spatial curvature, however, there is an additional source term in this 
equation, invalidating the argument that follows. 



(3.17) 




(3.18) 



= 2/3. 



(3.19) 
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This is only the case, however, for the Kaluza-Klein zero mode: all the higher modes 
have a rapid oscillatory time dependence such that the right-hand side does contribute at 
leading order. For these higher Kaluza-Klein modes, therefore, a gradient expansion does 
not exist. 

Setting (j) — 2(3 from now on, returning to (3.11) and recalling that (3 — 3 In 6, we 
immediately obtain the equivalent of the Priedmann equation on the branes, namely 

(t).+^(t)' = °- f^-^"* 

Integrating, the brane scale factors are given by 

6± = A±(r±-c±)^/^ (3.21) 

To fix the arbitrary constants A± and c±, we need only to expand the above in powers of 
T± and compare with (3.16). We find 

h±=[l± (4a tanh t/q) r±] ^'^ . (3.22) 

This equation determines the brane scale factors to all orders in in terms of the confor- 
mal time on each brane. (As a straightforward check, it is easy to confirm that the O(t^) 
terms in (3.16) are correctly reproduced). 

3.3 The scaling solution 

We are now in a position to solve the bulk equations of motion perturbatively in y^. 
Setting 

/3 = /3o + vl A + 0{yl) , v^vo + y>i + 0{yl) , (3.23) 

the leading terms of this expansion (namely (3q and uq) constitute a scaling solution, whose 
form is independent of yo for all y^ <^ 1. To determine the uj dependence of this scaling 
solution, we must solve the bulk equations of motion at zeroth order in y^. Evaluating 
the linear combination (3.6) — 3 [(3.7) + (3.8)], noting that at this order the right-hand 
sides all vanish automatically, wc find 

((/3^ - H^^y = 0. (3.24) 

Using the boundary condition P'^ — Suq = on the branes, we then obtain (3q = 3^0 + f{T), 
for some arbitrary function /(T). Substituting back into (3.7) and taking the square root 
then yields 

i/^ = ae-5'^«-2^ (3.25) 
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where consistency with the junction conditions (3.4) forced us to take the positive root. 
Integrating a second time, and re- writing e^^ = B~^^, we find 



^uo ^ qG^j^-j ^1/5^ g/3o 



(3.26) 



where h — bau + C{T), with C(T) arbitrary. 

In the special case where B and C are both constant, we recover the exact static 
domain wall solution (2.4), up to a trivial coordinate transformation. In general, however, 
these two moduh will be time-dependent. Inverting the relation h%. — S^(±5q; -|- C) to 
re-express B and C in terms of the brane scale factors h± (T) , we find 



B' = 



1 

10a 



6i), C = 5a 



bi-b^_ 



(3.27) 



Furthermore, at zeroth order in ^q, the conformal times on both branes are equal, since 



yodr=^dr = e--/^o/3^ = i3^(r)^ 

is independent of uj. To this order then, (3.22) reduces to 

6± = (l±4o;yoT)^/", 

allowing us to express the moduli in terms of r as 

1 



B'ir) 



10a 



C{t) = 5a 



[{1 + iayor f - {1 - iayorn 



{l+4ayorf/^ + {l~4ayorf/^ 



[1 + Aayor f/'^ - (1 - 4a|/o^)^/^ 
The brane conformal time r and the Milne time T are then related by 

InT = lOayo J [(1 + ^ayorf^ - (1 - AayorfT'^r. 



(3.28) 

(3.29) 

(3.30) 
(3.31) 

(3.32) 



Rather than attempting to evaluate this integral analytically, we will simply adopt r as 
our time coordinate^ in place of T. The complete scaling solution is then given by 

ds^ = h-'/\r,u;)\B\T){-dT^ + dx^) + B'\r)du;' 
e't' = S6(r)/i6/5(r,a;), 
h{T,uj) = 5acu + C(r), (3.33) 
''Note that the relation between T and t is monotonic. At small times T = y(jT-\-{a'^ /4)yQT^ +0{yQT^). 
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with B{t) and C(r) given in (3.30) and (3.31). This scahng solution solves the bulk 
Einstein equations and the junction conditions to leading order in our expansion in y^. 
The subleading corrections at higher order may be obtained in an analogous fashion, 
although we will not pursue them here. 

The expressions above for the scahng solution can also be used to calculate the slope 
of the 'warp factor' appearing in the metric (3.33), namely hB^. The slope, which from 
the Israel matching condition represents the effective strength of the two brane tensions, 
is given by 

ihB'Y = ^ [(1 + ^ayorf' - (1 - 4ayorf'] . (3.34) 

Thus we can picture the brane tensions to be evolving in time. In particular, note that 
what was a downward-pointing kink before the collision turns into an upward-pointing 
kink after the collision (and vice versa), while at the collision itself the slope is zero. 
Thus the tension of the branes swaps over at the collision, with both branes becoming 
effectively tcnsionless at the collision itself. (This is another indication that the collision 
represents a rather mild singularity). 

The scaling solution is valid for times in the range — (4ayo) ""^ < r < {4ayo)~^. At 
T = ±{Aayo)~^ , however, the scale factor on the negative-tension brane vanishes from 
(3.29) (recalling that the tension of the branes is reversed for r < 0). Since the physical 
interpretation of these events is much clearer in the alternative coordinate system in which 
the bulk is static and the branes are moving, we will postpone a full discussion until 
Section 4.1. (In fact, it will turn out that to continue the solution to times r > {iayo)'^, 
we simply need to introduce absolute value signs around all factors of (1 — Aayor)). 

A further quantity of interest, the distance between the branes, evolves as 

d = r h'/^'B'^ = ^ [(1 + 4a|/or)'/' - (1 - ^ayorf']. (3.35) 
7-1 oa 

Thus, for small r we have d ~ SyoT, independent of a, and for large r (imposing an 
absolute value on the second term), we have d ~ {yor/aY^'^. 

3.4 Lifting to eleven dimensions 

It is straightforward to lift the scaling solution to eleven dimensions, where the five- 
dimensional metric and scalar field are both part of the eleven-dimensional metric. The 
4-form field strength has a non-zero component in the Calabi-Yau directions only: 

Gabcd = —7= tabc/^ ujef, yo < y < yo, (3.36) 
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where cUg/ is the Kahler form on the Calabi-Yau. In the extended picture, the sign of 
Gabcd reverses across the branes. The metric in eleven dimensions reads 



= h-^/'B-^ (-dr^ + df 2) + h-'^^'B^du'' + 0'BWcy, (3.37) 

and the 4-form flux is proportional to a when all its indices are pointing in a Calabi-Yau 
direction. Thus the eleven-dimensional distance between the branes is given, to leading 
order in i/q, by 

r+l 

d^^ J dcu h-^/^ = 2yoT. (3.38) 

This simple relationship underlies the utility of r as a time coordinate. 

Finally, the volume of the Calabi-Yau manifold, averaged over the orbifold, is 



[I + AayoTfl^ - |1 - 4a|/or|5/4 



(3.39) 



where we have imposed the absolute values to permit a continuation to late times, as 
will be explained in Section 4.1. Thus the radius of the Calabi-Yau grows at large times 
as {e^l^) ~ T^/^, whereas the distance between the branes dn grows hnearly in r. A 
phenomenologically acceptable configuration where du ~ 30 e*^/^ [16] is therefore quite 
naturally obtained, assuming that both the distance between the branes and the Calabi- 
Yau volume modulus are stabilised by an inter-brane potential when they reach large 
values. We will leave a more detailed investigation of the eleven-dimensional properties 
of our solution to future work. 



4 An alternative perspective: moving branes in 
a static bulk 

Let us now consider an alternative coordinate system in which the bulk is static and the 
branes arc moving. To find the bulk metric in this coordinate system we can make use 
of a modified version of Birkhoff's theorem, as shown in Appendix A. Assuming only 
three-dimensional homogeneity and isotropy, as consistent with cosmological symmetry 
on the branes, in addition to the exact relation = 2/9 (motivated in Section 3.2 through 
considerations of regularity at the brane collision), the bulk Einstein equations can be 
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integrated exactly. One can then choose the static parameterisation: 

f{r) = a\^-nr\ e"^ = (4.1) 

where < r < a, and the coordinate t is unrelated to the Milne time appearing in the 
previous section. Physically, this solution describes a timelike naked singularity located 
at r = 0, and was first discovered by Chamblin and Reall in [10] (who instead looked 
directly for solutions in which the bulk was static). Since the coordinates above do not 
cover the whole spacetime manifold, the maximal extension may easily be constructed 
following the usual Kruskal procedure, as detailed in Appendix C. 

To find the trajectories of a pair of positive- and negative-tension branes embedded 
in this static bulk spacetime we solve the Israel matching conditions in the usual fashion. 
After performing this calculation (sec Appendix B, and also [10]), one finds that the 
induced brane metrics are indeed cosmological, with scale factors given by 

b^ = {l±A^r^f\ (4.2) 

where r± is the brane conformal time, and we have rescaled b± to unity at the collision, 
which is taken to occur at t± = 0. 

Upon setting fi = 0, we immediately recover the static domain wall solution (2.4), 
after a suitable change of coordinates. More generally, we require > to avoid the 
appearance of imaginary scale factors. Through comparison with our earlier result (3.22), 
we also find 

y/JI = a tanh^o- (4.3) 

Thus, for // > 0, after starting off coincident, the two branes proceed to separate. However, 
while the positive-tension brane travels out to large radii unchecked, the negative-tension 
brane reaches the naked singularity (at which 6_ and e'^ tend to zero) in a finite brane 
conformal time r_ = 

In the following section we will argue that the resulting singularity is extremely mild, 
and simple to regularise. If almost any type of well-behaved matter is present on the 
negative-tension brane - even in only vanishing quantities - then, rather than hitting the 
singularity, the negative-tension brane will instead undergo a bounce at some small finite 
value of the scale factor and move away from the singularity. 
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4.1 The bounce of the negative-tension brane 



The Priedmann equation for the negative-tension brane is derived in Appendix B. For 
the case in which a time-dependent scalar field (with an e*^ coupling to the Calabi-Yau 
volume scalar) is present on the brane, this equation takes the form (see (B.13)) 

^- = ^ = -2«-^ + -^ + ^, (4.4) 

where the constant A_ parameterises the scalar kinetic energy density. 

The key feature of this equation is the negative sign in front of the first term on the 
right-hand side, refiecting the fact that matter on the negative-tension brane couples to 
gravity with the wrong sign.^ For sufficiently large values of the scale factor the right- 
hand side is dominated by the jihZ^^ term. If we further assume that the matter density 
on the brancs is small compared to the brane tension^ {i.e. A_ ^ a), so that the term 
linear in A_ dominates over the quadratic term, then it follows that at some small value 
of the scale factor the entire right-hand side must vanish. Thus a negative-tension brane, 
initially travelling towards the singularity, will generically undergo a smooth bounce at 
some small value of the scale factor. After this bounce the brane travels away from the 
singularity back towards large values of the scale factor. This behaviour is specific to the 
negative-tension brane, and moreover, persists even in the limit in which A_ (and hence 
the initial matter density) is negligibly small. ^ 

In fact, even in the complete absence of matter on the negative-tension brane, the 
bounce off the naked singularity is still a relatively smooth event: converting to Kruskal 
coordinates so that light rays are at ±7r/4 angles, using our exact solution we show in 
Appendix C that the trajectory of the negative-tension brane becomes precisely tangential 
to the singularity at the moment of the bounce. In these coordinates then, the negative- 
tension brane merely grazes the singularity with zero velocity, before moving away again 
according to a well-defined smooth continuation. The complete solution is illustrated in 
Figures 1 and 4. 

'"'This property is a general feature of braneworld models, sec e.g. [17]. 

^As is in any case necessary for the existence of a four-dimensional effective description. 

''In the simplest version of heterotic M-theory, the scalars present on the negative-tension brane do 

not couple to (j) [7]. Nevertheless, in the limit of small matter density on the brane, it can be shown that 

the brane bounces in a manner identical to the case described above. When the scalars do not couple to 

(p, there are corrections to the bulk geometry in the vicinity of the bounce, but these become negligibly 

small as the scalar field density decreases to zero. 



16 



In this solution, the scale factor on the negative-tension brane is given exactly for all 
times post-collision by 

b^^\l-A^T^\'/\ (4.5) 

where r_ is the conformal time on the negative-tension brane. This result shows us how 
to continue our earlier scaling solution in brane-comoving coordinates (see Section 3.3) 
to times after the bounce of the negative-tension brane: we simply insert an absolute 
value sign into the expression (3.29) for 6_ (recalling that r is indeed the brane conformal 
time in these coordinates). The moduli B{t) and C{r) then follow from (3.27), which 
amounts to replacing all factors of (1 — Aayor) in (3.30) and (3.31) with |1 — 4q;7/ot|. It 
is easy to check that this continuation does not affect the smooth evolution of quanti- 
ties defined on the positive-tension brane, and that the range of the coordinate cu along 
the extra dimension is unaltered by this continuation (see Figure 2). Nevertheless, this 
brane-comoving coordinate system is not particularly well-adapted to deal with the actual 
moment of the bounce itself, as evidenced by the 'kinks' in Figure 2. With the presence of 
regulatory matter on the negative-tension brane, however, the bounce will occur at some 
small non-zero 6_, smoothing out these kinks. 

To check that the scahng solution is qualitatively correct near the bounce of the 
negative-tension brane, it is interesting to examine the leading behaviour of the exact 
bulk metric about the singularity. Setting = 5p + c, for some constant c, the metric 
(4.1) becomes 



ds^ = (5p + c 



,2/5 



-{a^ - /x(5p + c)2/5)dt2 + df2 + 



(4.6) 



— //(5p -|- c)^/^ 

which, after trivial re-scalings, to leading order in r reads 

ds^ = (5p + cf'[-dt^ + dx^ + dp% (4.7) 

This is precisely the form of our scaling solution (3.33) at a given instant in time. 

Note also that in the exact solution, the proper velocity of the negative-tension 
brane vanishes as it approaches the naked singularity: from (4.1) the proper veloc- 
ity is r^dr/(/dt), which scales as r'^dr/dt at small r. Yet from (C.5), dr/dt scales as 
f/r^ ~ 1/r^, and so the proper velocity tends to zero linearly with r. This result, that 
the proper velocity of the negative-tension brane tends to zero linearly with r as it ap- 
proaches the naked singularity, is also found in Kruskal coordinates, as shown in (C.14). 
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ayoT 




Figure 2: Continuation of the scaling solution (accurate to leading order in yo) beyond the 
bounce of the negative-tension brane at ay^T = 1/4, using the absolute value prescription. 
Contours of constant u are plotted against the fourth power of the three-dimensional scale 

factor b and the brane conformal time r. The contours shown are, from left to right, lo = — 1, 
—0.75, —0.5, 0, 0.5, and 1. Thus the left-most and right-most trajectories correspond to the 
negative- and positive-tension branes respectively (for which b\. = |1 zb AayQT\), with the region 
in between representing the bulk. 

5 Conclusions 

We have presented a cosmological solution describing the collision of the two flat boundary 
branes in heterotic M-theory. This solution is a significant step towards our goal of 
describing the cosmic singularity as a brane collision within the well-motivated framework 
of Hofava-Witten theory. Requiring the collision to be the 'least singular' possible, i.e., 
that the metric tends towards (2d compactified Milne)/Z2 x R3 times a finite- volume 
Calabi-Yau, has two important consequences. First, it selects a single solution to the 
equations of motion. Second, it shifts the singularity in the Calabi-Yau volume that one 
might have naively expected at the brane collision to two spacetime events before and 
after the brane collision. We have shown these two events to be very mild singularities, 
which are easily removed by including an arbitrarily small amount of matter (for example 
scalar field kinetic energy) on the negative-tension brane. Before the initial bounce of 
the negative-tension brane, and after the final bounce, the solution presented here can be 
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identified with that described by Chen et al. [11]. 

When the branes move at a small velocity, we expect to be able to accurately describe 
the solution using a four- dimensional effective theory (sec e.g. [18-22] and also [12]). We 
shall present such a description in a companion publication [23]. 

If our coUiding brane solution is to successfully describe the universe, we must also 
add potentials capable of stabilising the moduli; in particular the volume of the Calabi- 
Yau manifold, which determines the value of gauge couplings, and the distance between 
the branes, which determines Newton's constant of gravitation. These potentials also 
permit us to generate an interesting spectrum of cosmological perturbations. Although 
the required potentials cannot yet be derived from first principles, we can study the 
consequences of various simple assumed forms. The results will be presented elsewhere 
[24]. 
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A A modified Birkhoff theorem 

In this appendix we show how the assumption of cosmological symmetry on the branes, 
coupled with the exact result (p = 2(3 (derived in Section 3.2 following from our assumption 
that the Riemann curvature remains bounded as we approach the collision), is sufficient 
to uniquely determine the bulk metric. Moreover, with an appropriate choice of time- 
slicing, the bulk geometry can be written in static form. Our derivation closely parallels 
the analogous case for the Randall-Sundrum model [25] . 
We start with the fully-general metric ansatz 

ds^ = e^-^^/^{-dt^ + dy^) + e^^/W„ (A.l) 

where (j{t^y) and ilj(t,y) are arbitrary functions, and the coordinates t and y arc fresh, 
unrelated to those appearing earlier. To allow for the possibility of open, fiat or closed 
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brane geometries, we will take 

dnl = (1 - V)"' V + p'dnl, (A.2) 

where A; = — 1, Oor+1 respectively, and dQ"^ denotes the usual line element on a unit two- 
sphere. The above choice represents the most general bulk metric consistent with three- 
dimensional spatial homogeneity and isotropy, after we have made use of our freedom to 
write the {t, y) part of the metric in a conformally fiat form. 

In light of Section 3.2, we will additionally set (p = 2iIj. Then, evaluating the Einstein 
equations (3.5)-(3.9) in the bulk, after taking appropriate linear combinations we find 

'^,ty = (^,y'^,t + (^,t'^,y-^'^,t'^,y, (A.3) 

ttt + tyy = '^'^,t^,t + '2(^,yi',y-^^i-^i^l, (A.4) 

ttt-tw = -i'l + i^l + Qa\^''-^^^'\ (A.5) 

'y,tt-'y,yy = -f^ + f^ + a\''^-''^l' + 2ke^'^-'^l\ (A.6) 
Switching to the hght-cone coordinates, 

^=\{t-y), v = ]^{t + y), (A.7) 

these become 

^,«n = 2a„^,„ - 3^^„, (A.8) 

V',^^ = 2c7,,^ ,„ - 3V^J, (A.9) 

i^,uv = -V^,„V^,. + 6a2g2.-i4^/3^ ^A.IO) 

^,uv = -i^,utv + o?e'''-'''^" + 2ke"^-^^l\ (A.ll) 

Integrating the first two equations, we find 

e^-^^ ^i;^^V'{v)=i;,,U\u), (A.12) 

where U'{u) and V'{v) are arbitrary non-zero functions, and we will use primes to denote 
ordinary differentiation wherever the argument of a function is unique. It follows that 
and (7 take the form 

= ij{U{u) + V{v)), e^''-^^ = i^'U'V. (A.13) 

Consistency between (A. 10) and (A.ll) then requires A; = 0, and so only spatially flat 
three-geometries are permitted.^ Integrating once, we find 

^'e^ = 9/i - 9aV2'/'/3, (A.14) 



^Similar conclusions were reached by Chamblin and Reall in [10]. 
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where n is an arbitrary constant. 
The metric 

ds^ = -A,p'e''l'/MUdV + e2^/Mf^ (A.15) 
upon changing coordinates to 

r = e^/^ t = 3{V-U), (A. 16) 

then takes the static form 

(a^H — /xr*) 

The corresponding Calabi-Yau volume is given by 

e"^ = (A.18) 

B Brane trajectories 

In this appendix we embed a pair of moving branes into the static bulk spacetime (A. 17) 
and derive the Friedmann equations describing their trajectories. This is accomplished 
by solving the Israel matching conditions [15]: 

K^.--l{T^.-l9^,T^), (B.l) 

where and g^j^ denote the brane extrinsic curvature, stress tensor and induced 

metric respectively, and we are assuming a Z2 symmetry about each brane. 

Parameterising the trajectory of a given brane as t — T {tp) and r — R {tp), where tp 
is the brane proper time, the induced metric is 

ds^ = -dtl + R\tp) df^ (B.2) 

hence R may be associated with the relevant brane scale factor b±. The brane 4- velocity 
is then u°- — (T, R, 0), where, throughout this appendix, we will use dots to indicate 
differentiation with respect to tp. The constraint u"'Ua — —1 yields the additional relation 
T — {f + R^'^R^y/'^ / f . Similarly, the unit normal vector is given by 

na = ±{R'R, -J if + R''R'fl\ 0), (B.3) 

where the function /(r) is as defined in (4.1). (Note that the choice of sign corresponds to 
our choice of which side of the bulk we keep prior to imposing the Z2 symmetry. Keeping 
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the side for which r < R{tp) leads to the creation of a positive-tension brane, with normal 
pointing in the direction of decreasing r, requiring the positive sign for n^. Conversely, 
if we retain the r > R{tp) side of the bulk creating a negative-tension brane, the normal 
points in the direction of increasing r and we must take the negative sign for rio). 
The three-spatial components of the brane extrinsic curvature are then 

At = ^^nf = TR-\f + R''R'f" 9% (B.4) 

while the brane stress-energy is given by = =F6«e~'^5'^. The Israel matching condition 
then yields the Friedmann equation H^. — R^/R^ — h\/h\. — /i/b}^. Integrating, we find 
[10] 

b^^{l±5^t^y/', (B.5) 

where the constants of integration have been fixed by re-scaling b± to unity at the collision, 
which is taken to occur at — 0. In terms of brane conformal time t±, this relation 
reads 

b± = {l±A^T^y/\ (B.6) 

where the origins of r± are chosen so that the branes collide at r± = 0. 

This result may easily be generalised to scenarios in which matter is incorporated on 
one or both of the branes. For example, in the case where we add a scalar field x each 
brane, allowing for an arbitrary coupling to the Calabi-Yau volume modulus 0, the 
action should be augmented by the terms 

V^l{9xfF{4>). (B.7) 



/ 



Ud,u)=±l 

The brane stress-energy is now 

T± = TGae-'^g^, + ^^(0) [x,aX,b - \ 9^, [dx?]- (B.8) 

Since we are interested in cosmological solutions, we will consider the scalar field x to be 
a function of time only. Evaluating the three-spatial components of (B.l) then yields the 
modified Friedmann equation 



a_ x^ x^F'^ ji 
12 b% 242 + 



Hi = ±T7;^ + ^ + Tu^- (B.9) 



To find the dependence of x on the scale factor b±, it is necessary to evaluate the tptp com- 
ponent of (B.l). This yields the equivalent of the usual cosmological energy conservation 
equations. We start with 



Kt,t, - KabW'u' = u'^u'VaUb = -n.a", (B.IO) 
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where the acceleration a'^ = u'^V^u'^. Since a'^Uc = 0, the acceleration may also be written 
as a'^ — arf, where a = —K^tp- Then, since dt is a Killing vector of the background^, we 
have a — Ut/rit. Thus (B.l) reads 

Kft = ±^ -^{aR±- R'x'F) = ±^ - - x'F, (B.ll) 
tpt^ j^ej^ dip ^ 24 ^ ^ 24 ' ^ ^ 

where we have used (B.9). This leads to 

Rix'F),R = -12x'F, (B.12) 

i.e. we have x^F — 24A±i?~^^, for some constant A±. Thus the Priedmann equation 
becomes 

We must also satisfy the junction condition arising from the equation of motion. This 
is given by 

2 n«9„0± = Tl2a e^"^ - ^ x'i^,</>, (B.14) 

which leads to 

Consistency with (B.9) then requires the coupling to be 

F(0) = e"^. (B.16) 

It can be verified that this coupling is also consistent with the equation of motion for x- 
Another form of brane-bound matter that alters the Priedmann equations on the 
branes, while remaining consistent with the bulk geometry, is a cosmological constant A 
with e~'^ coupling to the Calabi-Yau volume scalar. In this case, the modified Priedmann 
equation is given by 

,-,-9 , « A 1 A^ it ,^ „, 

We will use the modified Priedmann equations (which are generally similar in form 
to (B.13)) in Section 4.1, to understand the bounce of the negative-tension brane in the 
vicinity of the singularity. 



^For a Killing vector an^ = aricS,'^ = acS,'^ = (_''u^\/bUc = u''\/b{ucS,'^) — u'^u^Si'b^c, where the last 
term vanishes by Killing's equation, V(6^c) = 0, hence anj = u''dbU^ = 
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C Kruskal extension of the bulk geometry 



In this appendix we consider the maximal extension of the bulk geometry (4.1), beyond 
the range < r < ^JJijoL for which / > 0. These calculations will lead us to the Kruskal 
diagram in Figure 1. Following Chen et al. [11], we start with the Eddington-Finkelstein 
coordinates u and v, defined via 

du^W{dt+jdr), dv^Aa^{dt-jdr). (C.l) 

The entire spacetime manifold is then covered by the Kruskal coordinates 

[/ = exp(^-^«^, V = exp(^^vy (C.2) 

where, for later convenience, we have introduced the constant A = ^JJijot = tanh|/o. In 
these coordinates, the metric then reads 

2 / (^UdV 2 1-2 Q^ 
Note that r should be understood here as a function of UV, given implicitly by 



UV = a exp 



2 / 1 3 1 , 



1 + f 




1 — f 


)1 



(C.4) 



where the re-scaled radial coordinate f = Ar, and a is defined such that cr = +1 for < 
f < 1, while (7 = — 1 for f > 1. (This choice generates the smooth continuation illustrated 
in Figure 3). Surfaces of constant f thus correspond to the hyperbolae UV — constant. In 
particular, we have chosen the constant of integration deriving from (C.l) such that the 
singularity at f = maps to the hyperbola UV = 1. The second constant of integration 
implicit in (C.l) may then be fixed by requiring that the time slice t = corresponds to 
the line U = V for < f < 1, in which case V/U = exp(2A^Q;^t). 

From the results of the preceding appendix, the trajectory of a brane (in the absence 
of any additional matter) is given, for < f < 1, by 

with the upper and lower signs corresponding respectively to a brane moving away from, 
and a brane moving towards, the singularity. Specialising to the case of two branes moving 
in opposite directions, upon integrating we find 

T2A^a't = f2 + ln(l - f') - A' - ln(l - A'), (C.6) 
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uv 




Figure 3: UV plotted as a function of f, according to (C.4) (taking A = 1). Note in particular 
the smooth continuation to f > 1 generated by the sign flip in a. The surfaces of constant f 

correspond to the hyperbolae UV = constant, with this constant being positive for < f < 1, 
and negative for f > 1. 



where we have fixed the constants of integration by requiring that the brane scale factors 
are unity at the coUision (hence f — \ aXt — Qi). 

Converting now to Kruskal coordinates, we find the trajectories of the branes are given 
parametrically, for all < f < oo, by 



U = T{f), 

u = g{f), 



V = Q{r) (outgoing); 

V — T{f) (infaUing) ; 



where the smooth functions 

1 



crexp 



exp 



(r + - + - r 



^ 3^1n|l-f|-^A2-iln(l-A2)) 



A3 

1 _ 1 f2 ^ 1 f3 _ ^^^^ + ^) + ^ ^2 ^ 1 ^^^^ _ y2^^ 



(C.7) 
(C.8) 

(C.9) 

(C.IO) 



As previously, (T = +lforO<f<l, but cr = — 1 for f > 1, ensuring that UV < for the 
portion of the trajectories parameterised by f > 1. 

To extend the brane trajectories beyond the bounce off the naked singularity, we first 
of all write down the 'extended' trajectories in terms of the original f and t coordinates. 
Altering the constants of integration in (C.6) so that the different branches of the trajec- 
tories match up at the bounce, we find 



T2A^a^t = + ln(l - f') + A^ + ln(l - A^), 



(C.ll) 
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for < r < 1, where the upper and lower signs correspond respectively to an outgoing, 
and an infalling, brane. Converting to Kruskal coordinates, this becomes 

U = T{f)g\Q), V = g{f)J^\Q) (outgoing); (C.12) 
U = g{f)T\Q), V = T{f)g\Q) (infalling); (C.13) 

where f now takes values in the entire range < r < oo. Finally, plotting (C.7), (C.8), 
(C.12) and (C.13), we obtain the full Kruskal diagrams shown in Figure 1 (for = 0.5), 
and in Figure 4 (for yo = 1 and 1.5). An interesting feature of the branc trajectories re- 
vealed by these plots are the points of inflection that occur whenever the brane trajectories 
intersect the Boyer axes UV = (see Figure 1 especially). These may be understood as 
a consequence of the vanishing of / in (C.5) whenever f = 1. 

Note also that at the singularity (f = corresponding to UV = 1), the slope of the 
negative-tension brane trajectory is exactly —l/U"^. For example, using (C.7), 

dl/ fl-f\V 

' ^ (C.14) 



dC/ \l + fj [/' 

which reduces to —1/U'^ when f = and UV = 1. A similar conclusion apphes for each of 
the other solution branches (C.8), (C.12), (C.13). Consequently, the brane trajectory at 
the bounce is tangent to the singularity itself. This means that the branc simply grazes 
the singularity with vanishing normal velocity. Similarly for the Calabi-Yau volume at 
the branes, Vcy = r^, we find 



dt 



±6A/ = ±6Aq;V + O(r^), (C.15) 



and hence the rate of change of the Calabi-Yau volume on the negative-tension brane, as 
it bounces off the singularity, is zero. The above considerations suggest that the bounce of 
the negative-tension brane is a relatively smooth event, even in the absence of regulatory 
matter on the brane. 
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Figure 4: Kruskal diagrams illustrating the brane trajectories for different collision rapidities. 
In the upper plot yo = 1, and in the lower plot i/q = 1.5. The corresponding plot for yo = 0.5 
appeared previously in Figure 1. 
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